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Abstract 

With the previous results for the analytical blowup solutions of the iV-dimensional (N > 2) 
Euler-Poisson equations, we extend the similar structure to construct an analytical family of 
solutions for the isothermal Navier-Stokes equations and pressureless Navier-Stokes equations 
with density- dependent viscosity. 



1 Introduction 



The Navier-Stokes equations can be formulated in the following form: 

Pt+V ■ (pu) =0, 

(1) 

(pu) t +V ■ (pu (g) u) + <5VP —vis(p, u). 

As usual, p = p(x, t) and u(x, t) are the density, the velocity respectively. P = P{p) is the pressure. 
We use a 7-law on the pressure, i.e. 

P(p) = Kp\ (2) 
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with K > 0, which is a universal hypothesis. The constant 7 = cp/c v > 1, where c p and c v are the 
specific heats per unit mass under constant pressure and constant volume respectively, is the ratio 
of the specific heats. 7 is the adiabatic exponent in ([2]). In particular, the fluid is called isothermal 
if 7 = 1. It can be used for constructing models with non-degenerate isothermal fluid. 5 can be the 
constant or 1. When S = 0, we call the system is pressureless; when S = 1, we call that it is with 
pressure. And vis(p,u) is the viscosity function. When vis(p 1 u) = 0, the system (fTJ) becomes the 
Euler equations. For the detailed study of the Euler and Navier-Stokes equations, see Q] and [I]. 
In the first part of this article, we study the solutions of the iV-dimensional (N > 1) isothermal 
equations in radial symmetry: 

Pt + Up r + pu r + ^-PU = 0, 

(3) 

p (ut + uu r ) + VKp =vis{p, u). 
Definition 1 (Blowup) We say a solution blows up if one of the following conditions is satisfied: 

(1) The solution becomes infinitely large at some point x and some finite time T; 

(2) The derivative of the solution becomes infinitely large at some point x and some finite time T . 

For the formation of singularity in the 3-dimensioanl case for the Euler equations, please refer 
the paper of Sideris [10]. In this article, we extend the results form the study of the (blowup) analyt- 
ical solutions in the TV-dimensional (N > 2) Euler-Poisson equations, which describes the evolution 
of the gaseous stars in astrophysics [2], [3], [7], [H] and [13], to the Navier-Stokes equations. For the 
similar kinds of blowup results in the non-isothermal case of the Euler or Navier-Stokes equations, 
please refer [S] and [T2J. 

Recently, Yuen's results in [13] , there exists a family of the blowup solution for the Euler-Poisson 
equations in the 2-dimensional radial symmetry case, 



The solutions are 



p t + up r + pu r + \pu = 0, 

p (u t + uu r ) + Kp r = JJ" p(t, s)sds. 



(4) 



fl ( < )= - sjjji «(°) = a o > °> a (°)= 01; 
y(x) + ly(x) + ^ e v(x) = p, y(0) = a, y(0) = 0, 



(5) 
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where K > 0, p, = 2A/K with a sufficiently small A and a are constants. 

(1) When A > 0, the solutions blow up in a finite time T; 

(2) When A = 0, if a\ < 0, the solutions blow up at t = — cto/cii- 

In this paper, we extend the above result to the isothermal Navier-Stokes equations in radial 
symmetry with the usual viscous function 

vis{p, u) = vAu, 

where v is a positive constant: 

p t + up r + pu r + ^-pu = 0, 

(6) 

p (U t + UU r ) + Kp r = v(u rr + ^-Ur - &^u), 

Theorem 2 For the N -dimensional isothermal Navier-Stokes equations in radial symmetry (0|), 
there exists a family of solutions, those are: 

p { t,r) = ^ey(r/«W, u (t,r) = $r, 

a(t) = o(0) = a > 0, a(0) = a u (7) 

y( x ) = m x2 + a > 

where a and A are arbitrary constants. 

In particular, for A > 0, the solutions blow up in finite time T. 

In the last part, the corresponding solutions to the pressureless Navier-Stokes equations with 
density-dependent viscosity is also studied. 



2 The Isothermal (7=1) Cases 



Before we present the proof of Theorem^ the Lemma 6 of [13] could be needed to further extended 
to the TV-dimensional space. 

Lemma 3 (The Extension of Lemma 6 of [13j) For the equation of conservation of mass in 
radial symmetry: 

N - 1 

Pt + Up r + pu r H pu — 0, (8) 

r 
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there exist solutions, 



a (tr ' [ ,r) ~a(ty 



P&r) = J -^A, «(t,r)=^r, (9) 



with the form f > e C 1 and a(t) >OeC 1 . 

Proof. Wc just plug |9]) into ©. Then 

N — 1 

Pi + Up r + pu r H pu 

r 

-Na(t)f(r/a(t)) a(t)rf(r/a(t)) 



ait)^ 1 a(t) N + 2 

a{t)rf(r/a(t)) f(r/a(t)) a(t) N - 1 / \r / a(t)) a(t) 
a[t) a(t) N + 1 a{t) N a(t) r a{t) N a(t) 

= 0. 

The proof is completed. ■ 

Besides, the Lemma 7 of [13] is also useful. For the better understanding of the lemma, the 
proof is given here. 

Lemma 4 (lemma 7 of |13j ) For the Emden equation, 



a(*) = -^ 
a(0) = ao > 0, d(0) = a\, 

we have, if A > 0, there exists a finite time T_ < +oo such that a(T_) = 0. 
Proof. By integrating (TIT)]) , we have 



(10) 



< ^a(i) 2 = -Alna(t) + 6> (11) 



where = Alnao + |af. 
From (fTTJl . we get, 

a(f) < e e/ \ 

If the statement is not true, we have 

< a(t) < e e/x , for all t > 0. 
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But since 



A -A 

a[t) = ~W)-^ 

we integrate this twice to deduce 

rt f T —\ — At 2 

a{t)<J J -^dsdr + dt + Co = ^jj- + dt + C Q . 

By taking t large enough, we get 

a(t) < 0. 

As a contradiction is met, the statement of the Lemma is true. ■ 

By extending the structure of the solutions (|5|) to the 2-dimensional isothermal Euler-Poisson 
equations (|4]) in [13] , it is a natural result to get the proof of the Theorem [2j 
Proof of Theorem [2j By using the Lemma[3j we can get that satisfy (j6])i. For the momentum 
equation, we have, 



By choosing 



p(ut + u ■ u r ) + Kp r — v(u rr H u r 

a(t) K r 
a(t) a{t) a[t) 

Pi Ar . r 
a(t) a(t) a(t) 



V(x) = ^x 2 +a, 



N-l N-l 



-u 



we have verified that ([7|) satisfies the above • If A > 0, by the Lemma IH there exists a finite 
time T for such that ct(T_) = 0. Thus, there exist blowup solutions in finite time T. The proof is 
completed. ■ 

With the assistance of the blowup rate results of the Euler-Poisson equations i.e. Theorem 3 
in [T3], it is trivial to have the following theorem: 

Theorem 5 With A > 0, the blowup rate of the solutions (7p is, 

( lim^,0)(T,-tr>O(l), 

where the blowup time T» and a < N are constants. 
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Remark 6 If we are interested in the mass of the solutions, the mass of the solutions can be 
calculated by:. 

/• r+oo 

M(t)= p(t,s)ds = a(N) p{t,s)s N - 1 ds, 
Jr n Jo 

where a(N) denotes some constant related to the unit ball in R N : a(l) = X; a{2) — 2ir; for N > 3, 



y(N) = N(N - 2)V(N) = N(N - 2)- 



r N/2 



T(N/2 + iy 

where V(N) is the volume of the unit ball in R N and T is the Gamma function. We observe that 
the mass of the initial time 0: 
(1) for A > 



M(0) 



a(N) r +0 ° 



a Q JO 



e 2K° + a S N - l ds. 



The mass is infinitive. The very large density comes from the ends of outside of the origin O. 
(2) for A < 0, 



, , a(N) r + °° 



2 o Jo 



+oc 



The mass of the solution can be arbitrarily small but without compact support if a is taken to be a 
very small negative number. 



Remark 7 Our results can be easily extended to the isothermal Euler/ Navier-Stokes equations 
with frictional damping term with the assistance of Lemma 7 in U2f : 



pt + up r + pu r + ^—^pu = 0, 



p(u t + u- u r ) + Kp r + f3pu = v(u rr + ^-^-Ur - ^pr-u) 



where [3 > and v > 0. 
The solutions are: 



P^ r ) = e "a'(t)^'" X*' r ) = 

a(t) + f3a(t) = ^4, a(0) = a > 0, a(Q) = a 1 



y{x) 
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2K 



xr + a. 



Remark 8 Our results can be easily extended to the isothermal Euler/ Navier-Stokes equations 
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with frictional damping term with the assistance of Lemma 7 in 

p t + Up r + pu r + ^-^-pu = 0, 

p(u t + u- u r ) + Kp r + (3pu = v(u rr + ^f^Ur - ^r-u) 
where [3 > and v > 0. 
The solutions are: 



p y(r-/a.{t)) 



a(t) + pa(t) = o(0) = a > 0, a(0) = oi, 

y( x ) = w^ 2 + a - 

Remark 9 The solutions {5p to the Euler-Poisson equations only work for the 2-dimensional case. 
But the solutions to the Navier-Stokes equations work for the N -dimensional [N > 1) case. 

Remark 10 We may extend the solutions to the 2-dimensional Euler/ Navier-Stokes equations 
with a solid core Jfjj/: 

p t + Up r + pu r + \pu = 0, 

p (u t + uu r ) + Kp r + (3pu _|_ v(u rr + \u r - ±;U), 

where Mq > 0, there is a unit stationary solid core locating [0, ro], where ro is a positive constant, 
surrounded by the distribution density. 
The corresponding solutions are: 

r. \ e H( '' /o(t)) n \ alt) t _ 

a(i) + /3d(i) = ^4, a(0) = a > 0, a(0) = a x , 
y{x) = -^x 2 + M Inx + a, 
where a > ^ is a constant. 



3 Pressureless Navier-Stokes Equations with Density-dependent 



Viscosity 



Now we consider the pressureless Navier-Stokes equations with density-dependent viscosity: 



vis{p,u) = v(m(p) V -u), 
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(12) 



in radial symmetry: 

Pt + up r + pu r + ^-f^-pu = 0, 

p (u t + uu r ) = (p(p)) r (^-u + u r ) + p(p)(u rr + - ^-u) 7 

where p(p) is a density-dependent viscosity function, which is usually written as p(p) = up with 
the constants k, 9 > 0. For the study of this kind of the above system, the readers may refer 



V 



(13) 



We can obtain the similar estimate about Lemma 2] to the following ODE, 

a(0) =a > 0, d(0) =ai < 
Lemma 11 For the ODE \13\) . with A > 0, there exists a finite time T- < +oo such that a(T_) = 

0. 



Proof. (1) If a(t) > and d(0) = a\ < for all time t, by integrating (fT3)l . we have 

A A A 

a W = 777 +0.1 < 777 • 

a(t) ao a(t) 



Take the integration for (TT 



t ,t 
a(s)a(s)ds < — / Ads, 



i[a(t)] 2 <-At + i a 2. 



When t is very large, we have 



~K*)] 2 <-i. 



A contradiction is met. The proof is completed. ■ 

Here we present another lemma before proceeding to the next theorem. 



(14) 



Lemma 12 For the ODE 



y(x)y(x) n -£x = Q, 
?/(0) = a > 0,n ^ -1, 



where £ and n are constants, 
we have the solution 



y (x) = - (n + l)& 2 + a n +K 



(15) 
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Proof. The above ODE (|15p may be solved by the separation method: 

y(x)y(x) n -£x = 0, 
y(x)y(x) n = £,x. 
By taking the integration with respect to x : 



y(x)y(x) n dx — / £xdx, 
ii Jo 



we have, 



y(x) n d[y(x)] = ^x 2 + C 1 , (16) 



where C\ is a constant. 

By integration by part, then the identity becomes 



y(x) n+ ' - n I y(x) n ~ 1 y(x)y(x)dx = ~£r 2 + d 



o 



1 - n I l/ (x)y(xrdx = ~Zx 2 + C 1 . 







From the equation (fT6|) . we can have the simple expression for y(x): 



.'/(•'• -n(^ !E a + C 1 ) = ^x 2 + r,. 



2/(x)" +1 = i(n+l)^ 2 + C 2 , 

where C2 = (n + 1)C\. 

By plugging into the initial condition for 2/(0), we have 

y (0) n+1 = =C 2 . 

Thus, the solution is: 



y{x) = n+ {j\{n+l)S,x 2 +a«+ 1 . 

The proof is completed. ■ 

The family of the solution to the pressureless Navier-Stokes equations with density-dependent 
viscosity: 



Pt + Up r + pu r + ^-r^-pu = 0, 



V 



(17) 

p (u t + uu r ) = (up ) r {^f 1 u + u r ) + np 9 (u rr + ^f^u r - ^r-u), 



is presented as the followings: 
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Theorem 13 For the pressureless Navier-Stokes equations with density-dependent viscosity Jt7| j 
in radial symmetry, there exists a family of solutions, 
for 6 = 1: 

p(t,r) = ^^-,u(t,r) - — • 



a W = a (°) = a o > 0, a(0) = ai, 



where a and A are arbitrary constants. 

In particular, for A > and a% < the solutions blow up in finite time; 
for 8^ 1: 

f ( 

°. Vy(^y)<o 

~ Aa(t) a(0) = a > 0, a(0) = a x , 




N9-N + 2 i 



where a > 0. 



Proof of Theorem 1131 To (fTT|) i . we may use Lemma[3]to check it. 
For 9 = 1, ifTTjlo becomes: 



(18) 



. . . ,N-1 . . N-l N-l . 

p[Ut + U ■ U r ) — (Kp)r( U + U r ) — Kp r (U rr H U r ^ U ) 

aft) Ke^M')) aft) 



aft) v a{t) N " a(t) 

/ Aa(t)r \ _ ^ (r/a(t)) y(^y) a(t) 
^a^J aft)- 

pa(t) / Ar 



Af+l 



(0 



aft) 2 Va(t) 



,2 i 777T -^( — ) 



>a(t) 



(19) 



where we use 



By choosing 



(HHJ) is equal to zero. 



Aot 
a(t) 2 



/ r s ■ / \ ^ 2 

y{~77\) — V\ x ) — 77~^~ x + a 



*(*) 
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For the case of 9 ^ 1, (fTT)) 9 can be calculated: 

. . g. (N-l \ g, N -I N-l . 

p[U t + U ■ U r j — [Kp ) r U + U r 1 — Kp (Urr H U r o U ) ' 



Xa{t)r \ NKOy^) -^^) a (t) 



a(t) m - N + 2 a(t) J a(t) N ( s ~ 1 )a(t) N + 1 a(t) 

' a(t) Ne - N + 2 a(t) ) a{t) N a{t) Ne ~ N + 2 

Xa{t )r \ N K ep y (^Y- 2 y(^)a(t) 

■ a (t) m - N + 2 a(t) J a(t) m ~ N + 2 



a(t) \ a(t) a(t) a{t) 



Define x n = 9 — 2, it follows: 

-pa(t) 



a(tY 



(Xx + NK9y(x) n y(x)^j (22) 



-AM*) ( x+ N*0 ymix)]t (23) 



and £ = -r^-* in Lemma [T2l and choose 



And this is easy to check that 

y(0) = 0. 

The equation (|2"2")) is equal to zero. The proof is completed. ■ 



Remark 14 By controlling the initial conditions in some solutions (18\j) , we may get the blowup 
solutions. And the modified solutions can be extended to the system in radial symmetry with 
frictional damping: 

p t + up r + pu r + ^-pu = 0, 



p (u t + uu r ) + (3pu = {p-{p))r (^-T~u + U r ) + p(p)(u rr + y^Ur - 
where (3 > 0, 

with the assistance of the ODE: 

a(0) = ao > 0, a(0) = ai, 
where S is a constant. 
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